A new type of vacancy- induced localized states in multilayer graphene 
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We demonstrate the existence of a new type of zero energy state associated to vacancies in 
multilayer graphene that has a finite amplitude over the layer with a vacancy and adjacent layers, 
and the peculiarity of being quasi-localized in the former and totally delocalized in the adjacent ones. 
In a bilayer, when a gap is induced in the system by applying a perpendicular electric field, these 
states become truly localized with a normalizable wavefunction. A transition from a localized to an 
extended state can be tuned by the external gate for experimentally accessible values of parameters. 

PACS numbers: 73.20.-r, 73.21.-b, 81.05. Uw 



Graphene is a one atom thick layer of carbon atoms 
ordered in a honeycomb lattice. The enormous interest 
risen since its discovery [l| is driven equally by potential 
technological applications [H and unconventional low- 
energy behavior (massless Dirac quasi- particles) To- 
gether with single layer graphene (SLG), bilayer graphene 
(BLG) and multilayer graphene (MLG) structures were 
also synthesized. The BLG structure being unique as a 
low energy effective model 0], rises better technological 
expectations due to the possibility to open and tune a gap 
in the spectrum by electric field effect [5|-|7| . A proper un- 
derstanding of the effect of disorder is crucial for technol- 
ogy relevant applications. Annealing and removing the 
substrate has recently lead to an increase in mobility by 
one order of magnitude in SLG [1, Intrinsic defects, 
such as vacancies or topological lattice defects are not 
easy to get rid of and further investigation on their role 
is mandatory. Vacancies are lately been recognized as 
one of the most important scattering centers in SLG and 
BLG [13, EH and the zero modes induced by this type 
of defects can greatly affect the transport properties of 
the samples as well as the possible electronic instabilities 
near the neutrality point. 

In the present manuscript we address the character of 
vacancy- induced electronic states in BLG, with extension 
to the case. We begin by summarizing the main findings 
of this work, (i) For the minimal model we construct an 
analytic solution on the lattice for the zero modes asso- 
ciated to the two different types of vacancies in BLG - 
located at A1/B2 or B1/A2 [see Fig.^a)]. We demon- 
strate that a new type of state different from these found 
in SLG and in other layered systems exists in BLG. The 
peculiarity consists in having a finite amplitude over the 
two layers and, more exotic, the wave function is quasi- 
localized in one layer and totally delocalized in the other. 
We also prove that these states survive in the continuum 
limit. The states associated to the Al/£?2-vacancies are 
quasi-localized decaying as 1/r at large distances, similar 
to these of SLG, while those corresponding to a B1/A2- 
vacancy are delocalized. The derealization is due to the 
spread of the wave function in the opposite layer where 
the vacancy resides. This solution is directly applica- 



Figure 1: (color online), (a) Bilayer lattice structure and 
main tight-binding parameters, (b)-(c) LDOS for 73 = and 
73 = Q.lt, respectively. 



ble to MLG and graphite, (ii) We demonstrate that 
these localization properties survive in the presence of 
non-minimal coupling 73 by means of a numerical anal- 
ysis including the study of the participation ratios, (iii) 
We study the behavior of these states in the presence 
of a gap and find that these associated to the B1/A2- 
vacancy become truly localized states leaving inside the 
gap while the AI/£?2-vacancy (monolayer type) become 
delocalized. The truly localized states are located sym- 
metrically around the middle of the gap - depending on 
the layer they belong to. 

Model. — The tight-binding minimal model for the 
7T— electrons in ^4i?-stacked BLG is schematically shown 
in Fig. QJa). We will use the parameters t ps 3 eV and 
71 ps t/10 [3j. It has two parabolic bands that touch at 
two degenerate Fermi points with a constant density of 
states (DOS) at the Fermi points. For the present study 
of zero modes it will also be of interest to consider the 
interlayer hopping 73 ps 71/3 that linearizes the bands 
around the Fermi points and induces a vanishing DOS at 
zero energy. Figure Q] shows the local DOS (LDOS) for 
73 = (b) and 73 7^ (c) . The presence of a finite gap in- 
duced through a perpendicular electric field E z = V/ (ed), 
where d ps 0.34 nm is the interlayer distance, is included 
by adding an on-site energy term: — V/2 at layer 1 and 
V/2 at layer 2. Within the present model vacancies cor- 
respond to the elimination of lattice sites. We do not 
include any reconstruction of the remaining structure. 
Even though some reconstruction might be present in 
real systems, the zero-energy modes we are interested in 
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Figure 2: (color online). Schematics for constructing a 
vacancy-induced zero-energy solution in bilayer graphene (see 
text). Circles indicate sites where the localized states have a 
finite amplitude, (a) A1/B2 vacancy, (b) B1/A2 vacancy. 



here seem to be rather insensitive to it |l2l |. 

Analytic construction of the vacancy states. — A va- 
cancy in the honeycomb lattice gives rise to a quasi- 
localized state fl3| whose wave function can be written 

as 



iK'. 



iy 



(1) 



where K and K' are the reciprocal space vectors of the 
two inequivalent corners of the first Brillouin zone, and 
(x, y) are distances in a reference frame centered at the 
vacancy position. We will construct an analytic solu- 
tion for vacancy-states in BLG within the minimal model 
following the analysis done for SLG in [l3|]. The wave 
function is obtained by matching surface state solutions 
at zigzag edges with those localized at Klein edges for 
a suitable boundary condition. The schematics used in 
this construction is shown in Fig. [5] The amplitude of the 
zero-mode wavefunction is denoted Ci(l,j), with c = a, b 
and i = 1,2 for sites in sublattice A, B and layer 1,2 
of the unit cell located at By cutting the system 

into left and right regions, defined with respect to the 
vacancy position, we see that in order to have a solution 
that decays away from the vacancy we need a zigzag- 
edge surface state to the left and a Klein-edge surface 
state to the right. The existence of surface states local- 
ized at zigzag edges in BLG has been proven in Ref. [3| ■ 



There are two linearly independent solutions, one living 
in a single layer (monolayer type) and the other having a 
finite amplitude over the two layers (bilayer type). Only 
sites belonging to the sublattice containing the zigzag 
edge have a finite amplitude. In an analogous way, one 
can show that surface states localized at Klein edges in 
BLG exist as well [l^. Again, two linearly independent 
solutions show up (monolayer and bilayer types). 

Consider a vacancy at A1/B2 sites, as sketched in 
Fig [2]^ a) . The zigzag- and Klein-edge surface states to 
be used have to have a finite amplitude, respectively, on 
the first zigzag column to the left and on the first beard 
column to the right of the vacancy. This is imposed by 
the matching (boundary) condition, which reads 



6 1 (-l,i) + 6i(0,i) + 6 1 (0,i-l) = 



(2) 



for all j's except at the vacancy, and involves sites of 
the two mentioned columns. The zigzag- and Klein- 
edge states with amplitudes starting at these columns 
are those of the monolayer type, i.e. with weight only on 
one layer: 

h(K0,3) = "£h(-l,k m )D-« + Ve ik ^+», (3) 



bi(l>0,j) = ^&i(0,M^l,e ifcm ' ( 



r+j) 



(4) 



where = — 2cos(fc/2), and the sums go over 27r/3 < 
k m < 4tt/3 in Eq. © and < k m , < 2tt/3 and 
47r/3 < k m i < 2tt in Eq. Q, for momenta k m , k m > along 
the y— direction, with bi(l,k) the Fourier transform of 
bi(l,j). The analysis now is completely analogous to the 
SLG case [^j. Namely, the boundary condition (|2]l. 
conveniently rewritten as J^k. 1, k m )e ^ = 

-Efc , (4 + elkm ' )M0, k m ,)e lk ™'^is satisfied for all k m 
and k m > in the ranges indicated above by choosing 
&i(— 1, k m ) = 1 and &i(0, k m >){l + e tk ™' ) = 1. Going from 
lattice indices (I, j) to distances (x,y) we obtain exactly 
the result given by Eq. ([1]). Therefore, for a vacancy at 
A1/B2 sites in BLG a quasi-localized (decaying as 1/r) 
zero-energy mode exists around the vacancy, living in the 
same layer but opposite sublattice. 

Consider now a vacancy at B1/A2 sites, sketched in 
Fig. Hfb). The zigzag- and Klein-edge states with a fi- 
nite amplitude, respectively, over sites (— and (0,j) 
of layer 1, are now those of the bilayer type. These 
states have amplitudes over layer 1 still given by Eqs. ([3]) 
and with the replacement b — > a. Additionally, they 
have also finite amplitudes over layer 2, which can be 
written as 



a 2 



«2 



(KOJ) = ^^ai(-4,fc m )a + l)^ +2) e ^ ( ^ + ^ 

(*>m - ^E«i(°' fc )^ + 1 )<' +1)et ' £m ' (i * i+J) '( 6 ) 
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with momenta k m ,k m i restricted to the intervals men- 
tioned before. An important point to note is that the 
boundary condition reads exactly the same as in Eq. ^ , 
with the replacement b — > a. Even though we are us- 
ing zigzag- and Klein-edge states which have finite am- 
plitudes in both layers, it happens that, by construction, 
the weight (O of the zigzag surface state at layer 2 is such 
that a2 (— = 0, and thus the matching condition at 
this layer is satisfied by default. At this point the deriva- 
tion follows closely that for a vacancy at A1/B2. Noting 
that in layer 1 we have to match exactly the same edge- 
state solutions given by Eqs. ([3]) and ((!]), with b — >• a, 
and that in layer 2 Eqs. ([5]) and © can also be written 
in the same form as Eqs. ([3J and apart from the term 
(/ + l)7i/t, we arrive at the following zero-mode behavior, 



T{x,y)~^{x,y) [1,171/*] 



(7) 



where ^(x, y) is the quasi- localized state given in Eq. ([lj, 
and the two component wave function refers to the two 
layers; first and second components for the first and 
second layers, respectively. This is a delocalized state, 
with the peculiarity of being quasi-localized in one layer 
(where the vacancy sits) and delocalized in the other 
where it goes to a constant when r — > oo. 

The analytic construction used for the minimal model 
in BLG applies directly to MLG and graphite with Bernal 
stacking along the lines of Ref. [TtJ ■ The quasi- localized 
state ((T|) is a solution in any multilayer with a A1/B2- 
vacancy. For a B 1/A2- vacancy the solution is a general- 
ization of state ([7]) with a quasi-localized component in 
the layer where the vacancy resides and delocalized com- 
ponents in the layers right on top and below this one: 
$(x,y) ~ V(x,y) [I,x7i/£,x7i/t]. 

The continuum limit. — Both the conventional [Eq. (fTJ)] 
and the unconventional [Eq. ([7])] solutions are fully 
consistent with the low-energy approximation for BLG 
[4[. Far from the vacancy the zero modes must obey 
dgipBi(z,z) = and d^ipA2{z,z) = at K, where 
z = x + iy and z = x — iy, and a similar set at K' 
with z replaced by z everywhere. An obvious solution 
has ipBi(z, z) = f(z) and ipA2(z, z) = 0, or 4>bi{.z, z) = 
and ipA2(z,z) = f(z), with f(z) analytic. Adding the 
contribution of the two if's we see that Eq. Q is pre- 
cisely of this form; the amplitude over the sublattice 
opposite to the vacancy behaving as 1/z + 1/z , analo- 
gous to the quasi- localized solution in SLG [13| . Inter- 
estingly, the bilayer model also supports solutions with 
ipBi(z,z) = zf(z) and i\)ai{z,z) = 0, or i/>bi(z,z) = 
and ipA2(z,z) = zf(z). Equation ([7J at the low-energy 
sublattice opposite to the vacancy is indeed a combina- 
tion of the stated solutions, namely z/z + z/z fl8l |. 

Vacancies in the gapless case. — The analytic results 
just presented are for 73 = 0. A finite 73 is crucial for 
the existence of the quasi- localized state (p}; otherwise a 
finite density of delocalized states exists in the same en- 
ergy region [see Fig. [IJb)-(c)]. This is addressed numeri- 
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Figure 3: (color online). LDOS (a)-(b) and IPR (c)-(d) for a 
vacancy at sublattice j41/_B2 (left panels) and B1/A2 (right 
panels). The LDOS is computed at a lattice site closest to 
the vacancy. The IPR is for the zero-energy mode induced by 
the vacancy. Lines are guides to the eyes. 



cally in the following. We also show that the delocalized 
character of the new solution (J7J persists in the presence 
of a finite 73. The localization character of vacancy- 
induced modes is studied through finite-size-scaling of 
the inverse participation ratio (IPR). The later is defined 
as V v = for the eigenstate v, where tp v (i) is 

its amplitude at site i. We perform exact diagonalization 
on small clusters with N up to 2 x 100 2 sites. The IPR for 
extended, quasi-localized, and truly localized states scales 
distinctively with N [16|. While for extended states we 
have V v ~ TV^ 1 , for quasi-localized states the 1/r decay 
implies V v ~ log(7V)~ 2 (consequence of the definition 
of the IPR in terms of normalized eigenstates) . For lo- 
calized wavefunctions the significant contribution to V v 
comes from the sites in which they lie, and a size inde- 
pendent V v shows up. Additionally to the IPR, we an- 
alyze the changes induced in the LDOS for sites around 
the vacancy. The LDOS is computed using the recursive 
Green's function method in clusters with ./V = 2 x 1400 2 , 
from which the thermodynamic limit can be inferred. 

In Fig. [3)Ja) and [31(b) we show the LDOS at a lattice 
site closest to a vacancy located in sublattice A1/B2 and 
Bl/^42, respectively. The sharp resonance at zero en- 
ergy in the former case is in agreement with the presence 
of a quasi-localized state, while the broader feature in 
the later may be attributed to the delocalized wavefunc- 
tion induced by a vacancy in B1/A2, which still presents 
a quasi-localized component in the layer where the va- 
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Figure 4: (color online). LDOS for a vacancy at sublattice 
Al (a) and Bl (b) for a finite gap, V = O.li. The LDOS is 
computed at a lattice site closest to the vacancy. Dashed lines 
are for the perfect lattice, (c) IPR averaged over the gap-edge 
resonance shown in (a), (d) IPR for the in-gap mode shown 
in (b). Dashed-dotted lines are guides to the eyes. 



cancy sits (and thus the feature). This interpretation is 
fully corroborated by the IPR scaling analysis shown in 
Fig. Etc) and Eld) for a vacancy in A1/B2 and B1/A2, 
respectively: quasi-localized state in the former case, and 
delocalized in the later. 

Vacancies in the gaped case. — When a finite electric 
field E z is present, a gap A g = [V 2 j 2 /(V 2 +^ 2 )] 1 ^ 2 opens 
between conduction and valence bands The quasi- 
localized state due to a vacancy at sublattice Al / B2 be- 
comes a resonance around ±V/2 in the gaped 
seen in the LDOS shown in Fig. Q|a) for a site closest 
to the vacancy. A strong resonance is seen around — V/2 
for a vacancy at Al (we used V = O.li), apart from the 
known gap edge divergence at — A ? /2 characteristic of 
the perfect lattice (dashed line) [19|. A vacancy at B2 
gives identical results with E — > —E. Such a vacancy- 
induced state living in the continuum is expected to be 
delocalized. This is confirmed by the IPR scaling as A -1 
(2fl | , as shown in Fig. |4jc) . 

For a vacancy at sublattice B1/A2, which originates 
the atypical delocalized state discussed above when no 
gap is present, a truly localized state inside the gap is 
induced when E z ^ 0. This is suggested by the sharp 
feature seen inside the gap in the LDOS for a site closest 
to the vacancy, as shown in Fig. EJb) and zoomed in the 
inset (marked by arrows). The IPR scaling to a constant, 
as seen in Fig.HJd), fully confirms the localized nature of 
this vacancy-induced state. Its asymmetric weight over 
the two layers explains why it appears off zero-energy: 



being negative for & Bl vacancy (layer 1 at an electro- 
static energy —V/2), as shown in Fig.[4jb), and positive, 
symmetrically placed with respect to the center of the 
gap, for a A2 vacancy (layer 2 at an electrostatic energy 
+V/2). 

Conclusions. — We have found a new type of zero mode 
state in BLG with special features: in the absence of a 
gap it is quasi-localized in one of the layers and delocal- 
ized in the other and in the presence of a gap becomes 
fully localized inside the gap. The results obtained in 
this work are directly applicable to MLG and graphite 
with Bernal stacking. The findings here reported can be 
important to understand recent experiments done in thin 
films of graphite irradiated with protons whose main ef- 



fect is to produce single vacancies on the sample [21 1. 
These samples show an enhanced local ferromagnetism 
that can be due to the local moments associated to the 
zero modes described in this work, and also a better con- 
ductivity than the untreated samples with less defects 
pointing to the idea that the delocalized states induced 
by the vacancies contribute to the conductivity. An en- 
hanced conductivity has also been found in acid-treated 
few-layer graphene (22|. The localized state found in the 
gaped case can also provide a natural explanation for 
the observation of localization inside the gap in the bi- 
ased BLG [?J and is in agreement with previous results 
obtained with impurity models in the continuum j23j. 
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